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Abstract: 

The results of a measurement of the proton structure function F2(x, Q 2 ) and the 
virtual photon-proton cross section are reported for momentum transfers squared 
Q 2 between 0.35 GeV 2 and 3.5 GeV 2 and for Bjorken-x values down to 6 • 10 -6 
using data collected by the HERA experiment HI in 1995. The data represent 
an increase in kinematic reach to lower x and Q 2 values of about a factor of 
5 compared to previous HI measurements. Including measurements from fixed 
target experiments the rise of F2 with decreasing x is found to be less steep for the 
lowest Q 2 values measured. Phenomenological models at low Q 2 are compared 
with the data. 
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1 Introduction 



The measurement of the inclusive deep inelastic lepton-proton scattering (DIS) cross section 
has been of great importance for the understanding of quark and gluon substructure of the 
proton O. Recently published measurements of the structure function F% at HERA based 
on analyses of the data collected in 1994 cover the range in squared four-momentum transfer 
Q 2 , Bjorken-x and inelasticity y corresponding to 1.5 < Q 2 < 5000 GeV 2 , 3- 10 -5 < x < 0.32 
and roughly y > 0.01 |2|-44|. These data have shown that i*2 continues to rise strongly with 
decreasing x for the lowest Q 2 values reached. Furthermore, it has turned out that the data 
can be successfully described by perturbative QCD in the measured kinematic range, using 
the leading twist Next to Leading Order QCD evolution equations, which have subsequently 
been used to extract the gluon density in the proton at low x. 



The measured behaviour of F2 at low Q 2 is contrary to expectations based on Regge 
phenomenology ||, which anticipate a much slower rise with decreasing x. At low x the 
center of mass energy W of the 7*p system is approximately W — \/Q 2 /x. Since F2 at low x 
is directly proportional to the total cross section of the virtual photon-proton interaction 
the rise of F^ with decreasing x is reflected in a strong increase of <t*°* with W. The rise can 
be quantified by parameterizing the data in the form F% ~ W 2X . For low-x data with Q 2 > 1.5 
GeV 2 and 40 < W < 200 GeV, values of A in the range of 0.2-0.4 have been measured [01. 
This can be contrasted with the increase of the real photoproduction (Q 2 = 0) cross section 
in the same range of W |], [?J which gives a value of A ~ 0.08, similar to that observed 
for hadronic interactions. These processes are dominantly of non-perturbative nature and 
the total cross section behaviour with W is well described by Regge phenomenology which 
assumes that the interaction dynamics is driven by the so-called soft pomeron. The F2 data 
measured in the region Q 2 > 1.5 GeV 2 indicate that the photoproduction limit has not 
yet been reached. It is therefore of interest to study data at still lower Q 2 values, and to 
experimentally establish and probe the transition between the region of perturbative QCD 
(DIS) and Regge phenomenology (photoproduction), a topic of much theoretical debate at 
present (see e.g. 



The data presented in this paper extend the Q 2 region in a continuous way down to lower 
values. This extension was achieved with the HI detector at HERA by upgrading the detector 
components in the backward (electron]] beam direction) region. A new calorimeter and a new 
tracking chamber with increased acceptance for small scattering angles were installed in the 
winter shutdown 1994/1995, and they were commissioned during the data taking period of 
1995. In 1995 the incident electron energy was E e = 27.5 GeV and the proton energy was 
E p = 820 GeV, leading to a total center of mass energy of the collision of y/s = 300 GeV. The 
data used for this analysis result from a short data taking period where the ep collision vertex 
was shifted by 70 cm in the proton-beam direction with respect to the nominal vertex position 
to increase the acceptance for the electrons scattered through small angles. As a result the 
data have an acceptance in Q 2 down to Q 2 ~ 0.3 GeV 2 . In this paper we present data for the 
kinematic range 0.35 < Q 2 < 3.5 GeV 2 , x > 6 • 10~ 6 and 0.03 <y< 0.75, and compare the 
results with other experiments and with Regge and QCD inspired phenomenological models. 

1 HERA operated with e~p collisions in 1992, 1993 and the start of 1994, and e + p collisions for the major 
part of 1994 and all of 1995. In this paper the incident and scattered lepton are always referred to as "electron" . 
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2 The HI Detector 



The HI detector p| is a nearly hermetic multi-purpose apparatus built to investigate the 



inelastic high-energy interactions of electrons and protons at HERA. The structure function 
measurement relies essentially on the tracking chamber system, the backward calorimeter 
and the liquid argon (LAr) calorimeter which are described briefly here. 

The tracking system includes the central tracking chambers, the forward tracker mod- 
ules and a backward drift chamber. These chambers are placed around the beam pipe at z 
positions between -1.5 and 2.5 m. The +z-axis is in the proton beam direction. A supercon- 
ducting solenoid surrounding both the tracking system and the (LAr) calorimeter provides a 
uniform magnetic field of 1.15 T. 

The central jet chamber (CJC) consists of two concentric drift chambers covering a polar 
angle range from 15° to 165°. Polar angles are defined with respect to the +z direction. Tracks 
crossing the CJC are measured with a transverse momentum resolution of 5pt/pt < 0.01 -px 
(GeV). The CJC is supplemented by two cylindrical drift chambers at radii of 18 and 47 cm, 
respectively, to improve the determination of the z coordinate of the tracks. A proportional 
chamber is attached to each of the z drift chambers for triggering. 

A tracking chamber system made of three identical modules measures charged particles 
emitted in the forward direction (7° to 20°). The forward tracker is used to determine the 
event vertex for the events which leave no track in the CJC. This enables vertex reconstruction 
of events with larger x than can be achieved with the CJC alone. 

In the backward region there is an eight layer drift chamber (BDC) which has a 
polar angle acceptance of 151° to 177.5° for collisions at the nominal vertex position. It 
replaces the backward multiwire proportional chamber used in previous data taking periods. 
For data where the vertex is shifted by 70 cm in the proton direction the acceptance at 
large angles increases to a maximum value of 178.3°. The BDC provides track segments for 
charged particles entering the backward calorimeter. These are used to identify electrons and, 
together with the event vertex reconstructed from tracks in the forward and central tracker, 
to measure their polar angle 6 e . The resolution for reconstructed BDC hits is about 0.5 mm 
in the radial direction, and 2.5 mm in the azimuthal direction. 

In the backward region a lead/scintillating fiber calorimeter (SPACAL) |]l4| was in- 
stalled, replacing the previous lead/scintillator electromagnetic calorimeter (BEMC). The 
new calorimeter has both electromagnetic and hadronic sections. The angular acceptance 
of the SPACAL is 153° < 6 < 177.8° for collisions at the nominal vertex, increasing to 
178.5° for collisions where the interaction vertex is shifted by 70 cm in the proton direc- 
tion. The resolution in electron energy is determined using the present data and a value of 
7.5%/-\/ -E(GeV) © 2.5% is obtained for the electromagnetic calorimeter. The absolute energy 
scale uncertainty is determined to be 1% at 27.5 GeV increasing to 3% at 7 GeV, which is 
the lowest electron energy used in this analysis. The high granularity (1192 cells) results 
in a spatial resolution of about 4 mm. Around the beam pipe the electromagnetic section 
of the SPACAL has a veto layer (inner and outer radius 5.7 cm and 6.5 cm respectively), 
which is used to detect electrons for which the shower is not fully contained laterally. The 
hadronic section has 128 cells. The hadronic energy scale uncertainty of the measurement 
in the SPACAL is presently about 7%. Both calorimeter parts have a time resolution better 
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than 1 ns which facilitates the reduction of proton beam induced background from beam-wall 
and beam-gas interactions which occur upstream of the detector. 

To summarize, compared with the detector used up to the end of 1994, the upgrade of 
the detector components in the backward region has resulted in an increased acceptance of 
electrons scattered through small angles (large 9 e ), in improved granularity and in improved 
resolution. When taken together with the hadronic energy measurement these features have 
greatly facilitated the recognition and removal of background when selecting DIS events, and 
have extended the available kinematic range of the measurement. 

Hadronic final state energies are further measured in the LAr calorimeter which covers 
an angular region between 4° and 154°. The calorimeter consists of an electromagnetic section 
with lead absorber plates and a hadronic section with stainless steel absorber plates. Both 
sections are highly segmented in the transverse and longitudinal directions with about 44000 
cells in total. The electromagnetic part has a depth of between 20 and 30 radiation lengths. 
The total depth of both calorimeters varies between 4.5 and 8 interaction lengths. The 
hadronic energy uncertainty of the LAr calorimeter is determined to be 4%. 

The luminosity is determined from the measured cross section of the Bethe Heitler (BH) 
reaction ep — ► ep^f. The final state electron and photon are detected in calorimeters (electron 
and photon "taggers" ) close to the beam pipe but at large distances from the main detector 
(at z = —33 m and z = —103 m). For the final value of the luminosity only the hard photon 
bremsstrahlung data are used. The precision of the luminosity determination for these data, 
where the interaction vertex is shifted, is 3% [16|. It results from the error on the luminosity 
measurement and a correction for proton "satellite" bunches in HERA which lead to collisions 
about 70 cm displaced with respect to the main part of the proton bunch. 



3 Kinematics 

The kinematic variables of the inclusive scattering process ep — ► eX can be reconstructed in 
different ways using measured quantities from the hadronic final state and from the scattered 
electron. The choice of the reconstruction method for Q 2 and y determines the size of 
systematic errors, acceptance and radiative corrections. The methods used in the analysis of 
the 1995 data are: i) the "electron method" (E), which uses only the event vertex and the 
reconstructed scattered electron, and which has the best resolutions in x and Q 2 at large y; 
ii) the "S method" (E) [17| , which uses the electron and hadronic final states measurements, 
which is less sensitive to radiative corrections, and which can be used from very low to large 
y values. Both calorimetric and track information are used to calculate the kinematics with 
the T, method [Q. The application of different methods is an important cross check of the 
results. 

The basic formulae for Q 2 and y in the E method are: 

1 E 'e ■ 2 #e ^2 a T?' TP 2 ® e E 2 sm 2 e 

y e = l- -sin - Q e = AE e E e cos — = — , (1) 

bj e l l l — y e 

where E' e and 9 e are the energy and polar angle of the scattered electron. The formulae for 
the X method are: 

„ - S n 2 _ E /2 sm 2 e e 
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with 

E = ^(E h -P Z;h ). (3) 
h 

Here Eh and P z ,h are the energy and longitudinal momentum component of a particle h, 
the summation is over all hadronic final state particles, and masses are neglected. The 
denominator of is equal to twice the energy of the true incident beam energy, which differs 
from the nominal beam energy if one or more real photons are (mostly collinearly) emitted 
by the incident electron and not detected. The variable x is calculated in both methods as 
x = Q 2 /ys. The methods used are the same as those in the analysis of the 1994 data which 
is reported in [||. 



4 Monte Carlo Programs and their Implementation 

Acceptance corrections and background contributions are studied with the data and Monte 
Carlo simulations. Monte Carlo DIS events corresponding to twice the luminosity of the 



data were generated using DJANGO [18]. This program is based on HERACLES [19| for the 



electroweak interaction and on LEPTO [20| to simulate the hadronic final state. HERACLES 
includes complete first order radiative corrections, the simulation of real bremsstrahlung 
photons, and the longitudinal structure function. LEPTO uses the colour dipole model 



(CDM) as implemented in ARIADNE [21] which is in good agreement with data on the 



energy flow and other characteristics of the final state as measured by HI 22] and ZEUS [23]. 



Alternatively, first order QCD matrix elements with additional parton showers can be used 



for the final state QCD radiation. Hadronization is performed using string fragmentation [24]. 
This model does not contain events with large rapidity gaps J25|, tentatively interpreted as 
diffractive events. Such events can be generated with the Monte Carlo programs RAPGAP 



[26] and DIFFVM |p7| ]. The latter generates the diffractive exclusive channels ep — > epp and 



ep — > ep<p. In the present analysis these channels to a large extent escape selection. Dedicated 



measurements of vector meson production cross sections have been made at HERA [28,29], 
and are used to normalize the Monte Carlo prediction. The program RAPGAP generates 
events with a continuous mass spectrum of diffractive final states and the yield has been 
normalized to the rate of events with a large rapidity gap observed in the data. 



The acceptance corrections were performed using in turn the GRV [30 and MRSDO' [32] 
parton distributions for the initial structure functions in the Monte Carlo calculations. An 
iterative procedure was used to reweight the input structure functions of the Monte Carlo 
programs with the measured F2 values in this analysis, as described in Section 6. For the 
figures in this paper which show a comparison of the detector response of Monte Carlo with 
data, the result of the final iteration of the input structure function was used. 

Photoproduction background was simulated using the PHOJET |33| generator for jp 
interactions. A sample of photoproduction events was generated which contained all classes 
of events (soft hadronic collisions, hard scattering processes and heavy flavour production), 
corresponding to three times the luminosity of the data. PHOJET was used to generate 
events with Q 2 < 0.1 GeV 2 , and DJANGO for events with Q 2 > 0.1 GeV 2 for acceptance 
and background calculations. The results reported here were found to be insensitive to the 
value of the Q 2 boundary between the regions where the two generators were used. 
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The detector response of the Monte Carlo events was simulated in detail with a program 
based on the GEANT program [34]. After this step these events were subjected to the same 
reconstruction and analysis chain as the real data. 



5 Event Selection and Calibration 

The data sample used for this analysis corresponds to an integrated luminosity of 114 nb _1 . 
The trigger used requires that there be a local energy deposit (cluster) in the SPACAL 
calorimeter with energy greater than about 5 GeV occurring in time with an ep bunch cross- 
ing. The trigger efficiency is about 99% for electrons with an energy above 7 GeV. Losses of 
about 1% occur due to the event timing requirements. 

Deep inelastic events are selected if they satisfy the following criteria: 

• the most energetic cluster in the electromagnetic section of the SPACAL is an electron 
candidate (see below) with a signal within a time window of 10 ns total width around 
the expected value for a genuine ep collision; 

• a reconstructed event vertex exists, as determined using the central or forward trackers, 
within 30 cm of the average event vertex position along the beam (z) direction; 

• for the electron method — P Z h) is required to be larger than 35 GeV, where Eh 
and P z h are the energy and longitudinal momentum of a particle. The sum is over all 
energy deposits measured with the calorimeters. 

For the electron candidate the following is required: 

• the energy of the cluster must be larger than 7 GeV; 

• the radius of the cluster^ must be smaller than 3.5 cm; 

• there must be a track segment in the BDC matched to the cluster in the SPACAL 
within 2 cm along the radial direction and within 2.5 cm along the azimuthal direction; 
the segment closest to the center of the cluster is used for the calculation of 9 e ; 

• the radial distance from the beamline to the point at which the track associated with the 
cluster intersects the surface of the SPACAL must be larger than 8.7 cm, corresponding 
roughly to an acceptance of 9 e < 178°; 

• less than 1 GeV of energy should be deposited in the veto layer of the SPACAL to avoid 
having too much loss of energy near the beam pipe; 

• the energy measured behind the electron cluster in the hadronic part of the SPACAL, 
within a radius of 17.5 cm of the projected electromagnetic shower center, should be 
less than 0.5 GeV. 

2 The cluster radius is defined as - Ei * di/ ^ . Ei, where the sum is over all cluster cells; Ei is the energy 
of cell i, and d 4 is the distance from the cluster center of gravity to the center of the cell i. 
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The selection cuts were designed to have a high efficiency for detecting DIS events. For a 
large part of the kinematic region studied the total efficiency is better than 90%. 

The main non-ep backgrounds in the event sample selected are due to proton beam 
interactions with residual gas and beam line elements upstream of the HI detector. An 
efficient reduction of the background is provided by the minimum electron energy and the 
vertex requirements discussed above. The residual non-ep background was estimated by 
visual inspection to be less than 2% of the total number of events at the highest y, and 
negligible elsewhere. 

The only significant background to DIS from ep interactions is due to photoproduction 
events (Q 2 ~ 0) in which the scattered electron escapes the main detector along the beam 
pipe and the electron signal is faked by an energy deposition associated with the hadronic final 
state. For about 10% of these events the scattered electron is detected in the electron tagger, 
and such events can be identified as photoproduction background. The total photoproduction 
background was estimated from studies using simulated events from the Monte Carlo program 
PHOJET. The results are presented in Fig. [I], where the energy spectrum of the downstream 
electron tagger and the energy of the (fake) electron candidate in the SPACAL are shown for 
tagged photoproduction events in the data and in the Monte Carlo predictions. The Monte 
Carlo curves are normalized to the number of tagged events. The prediction for the remaining 
photoproduction background was subtracted statistically bin by bin. For each Q 2 value, only 
the lowest x bin has a contamination larger than 5%. This contamination never exceeds 20%. 

The energy scale for electrons has been determined with events at low y, for which the 
energy of the scattered electron is very close to the incident electron energy. The linearity of 
the energy response was verified with QED-Compton events for the energy range used in this 
analysis. The precision of the angular measurement from the BDC was estimated using tracks 
in the central tracker extrapolated into the BDC region. Fig. ^ shows the energy and angular 
distribution of electron candidates in the sample of selected DIS events, together with Monte 
Carlo predictions for the sum of DIS and photoproduction background in the kinematic range 
0.32 < Q 2 < 10 GeV 2 . The structure function of the Monte Carlo calculation is reweighted 
as described in Section 6. For these and following figures in this section the Monte Carlo 
prediction is normalized to the luminosity. Good agreement is observed between data and 
Monte Carlo calculations. Fig. || shows the distribution of £(£tj — P z ,h) f° r the full data 
sample and for the sample with y > 0.55. The Monte Carlo calculation describes the data 
well. The enhancement observed around 25 GeV in Fig. |||b results from events with one or 
more photons collinearly emitted by the incident electron. The level to which the hadronic 
variables are understood is demonstrated in Fig. ||, where the distributions of the ratio of 
the y values measured with the E and E methods, yT,/y e , and the transverse momentum 
of the hadronic system and the scattered electron, Pt,h/pt,e, are shown. Here the data are 
limited to > 0.05 to ensure good quality reconstruction of the kinematics with the electron 
method, and to the region where the X method is used for the combined measurement, namely 
0.75 < Q 2 < 4.2 GeV 2 . In all, the detector response is well understood, allowing a precise 
extraction of the cross section and F2. 
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6 Structure Function and Cross Section Measurement 



The measured ep cross section in the HERA kinematic range can be expressed in terms of 
proton structure functions or cross sections for virtual photon-proton interactions as follows 

d 2 a 2-na 2 , y 2 . „ . ~ N 



dxdQ 2 Q 4 x" u l + R' zv ' (4) 

r ef J( 



T[a T (x, Q 2 ) + e(y)a L (x, Q 2 )} = Fat/Mx, y, Q 2 ). 



Here R = Fl/ {F^—Fl) where Fl is the longitudinal structure function, a is the fine structure 
constant, and ctl and ctt are the cross sections for transverse and longitudinally polarized 
virtual photons. The flux factor, V, and the ratio of the longitudinal to the transverse flux, 
e, are taken to be 

a{2-2y + y 2 ) 2(1 - y) 

r = 2nQ 2 x ' e{v)= 2-2y + y 2 - (5) 

The quantity a^lp is the effective measured virtual photon-proton cross section for ep colli- 
sions in our kinematic range, and can be determined from the data without assumptions for 
R. The total virtual photon-proton cross section is here defmedP] as 



47r 2 a 

~Q 3 



a% = a T (x, Q 2 ) + a L (x, Q 2 ) ~ ^-F 2 (x, Q 2 ). (6) 



With this definition cr*°' p depends only on Q 2 and x (or W) and the results of different 
experiments may easily be compared. 

The virtual photon-proton cross section is determined by converting the measured number 
of events in a given bin into a bin averaged cross section using Monte Carlo acceptance 
calculations. The data are binned in a grid of Q 2 and x for the region Q 2 > 0.75 GeV 2 
as for previous HI analyses, and a grid of Q 2 and y for the region below 0.75 GeV 2 , which 
optimizes the access to the smallest possible (Q 2 ,x) values. Bin widths are chosen such that 
the number of events reconstructed in any given bin which originate from that bin is larger 
than 40% for the E method and larger than 30% for the £ method. All detector efficiencies 
are determined from the data utilizing the redundancy of the apparatus. Apart from small 
extra corrections, all efficiencies are correctly reproduced by the Monte Carlo simulation, and 
therefore the Monte Carlo can safely be used to correct for acceptance and efficiency effects. 
The bin averaged cross section is corrected for first order QED radiative contributions with 
the program HERACLES. The effective virtual photon-proton cross section, <J^lp, is finally 
obtained by correcting the bin averaged cross sections for each bin to the values at the given 
bin centers. 

To extract the structure function F% from these measurements an assumption has to be 
made for the longitudinal structure function Fl since it has not yet been measured in this 
kinematic region. In this analysis the model of |37j is applied for the calculation of R. These 
values of R are then used to determine the F% values, as well as to reweight Monte Carlo 
events. The model is based on the photon-gluon fusion process and has the proper limit 



3 The exact formula used is a*J^ p = (47r 2 a/Q 4 )(4MV + Q 2 )/(l - x) ■ F 2 (x, Q 2 ) [§§§§, which is approxi- 
mately equal to eqn. ^ in the HERA kinematic range; M is the proton mass. 
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for Q 2 — ► where Fl should vanish oc Q 4 . The predictions of this model for the values 
of R in our kinematic region vary from 0.1 at Q 2 = 0.35 GeV 2 to 0.3 at Q 2 = 3.5 GeV 2 . 
The predictions for R at higher Q 2 are in agreement with measurements from fixed target 
experiments. It should be noted that this is a model and future measurements could reveal 
quite different R values. However, only the lowest x point at each Q 2 is affected significantly 
by the assumption made for R. If R is taken to be zero rather than the values obtained using 
the above model, the variation in F 2 is 5 to 10% at the highest y (smallest x) at a given Q 2 , 
and smaller elsewhere. 

Values of F2 are derived using an iterative procedure. An initial determination of the F2 
values uses a structure function parameterization in the Monte Carlo simulation as described 
in Section 3 and it assumes a value for R as given above. A new structure function for 
the full range in Q 2 is then calculated following the prescription of Badelek and Kwiecinski 
(BK, |J38| ), where the structure function is assumed to be the sum of two contributions: a 
Vector Meson Dominance (VMD) model term F^ MD a nd a partonic term F 2 art . The latter 
becomes dominant above Q 2 ~ 1 GeV 2 . In this analysis the result of a QCD fit similar to 
that reported in || is used for the partonic term. This fit was made to structure function 
data from HI M, and in order to constrain the high-x region, to NMC and BCDMS ]40| 
data. Parton density parameterizations were defined at a starting scale Qq = 0.35 GeV 2 , 
and data with Q 2 > 1 GeV 2 were fitted, yielding values for F2 denoted as p^iQCD ^ ^ e 
following. The newly measured low Q 2 data points and the measurements at Q 2 = in the 
W range at HERA |,| were fitted to the form 

F 2 (x, Q 2 ) = C VM F 2 VMD (x, Q 2 ) + n2 f n2 F 2 mQCD (x, Q 2 + Qv M ), (7) 

with x = (Q 2 + Q\/m)/(^ 2 + Q 2 + Qvm)- The fit parameters are the meson mass cut-off 
parameter Qy M and the normalization of the vector meson term Cvm- The latter parameter 
is not part of the BK model and was introduced to reproduce the real photoproduction data 
measured at HERA in a phenomenological way. Values of Q V m = 0.45 GeV 2 and Cvm = 0.77 
are obtained^. It was checked that no further iteration step was needed. The Monte Carlo 
curves discussed in Section 4 are reweighted with these F 2 values (and the R values discussed 
above). 

A list of sources of systematic errors in the F 2 determination is given below. 



• Uncertainty of the electron energy scale in the SPACAL, varying from 1% at large 
electron energies to 3% at 7 GeV. 

• A 4% scale error for the hadronic energy in the LAr calorimeter, the effect of which is 
reduced due to the joint consideration of tracks and calorimeter cells for the S analysis. 
A 7% scale error was assigned to the energy of the hadronic final state measured in the 
SPACAL. 

• A potential shift of up to 0.5 mrad for the electron polar angle. 

• For the electron identification efficiency the error was taken to be 30% of the fraction 
of events lost by the cuts, as given by the DIS Monte Carlo. 

4 Note that this procedure does not guarantee a consistent separation of F2 into the p^ MU an j pp art 
contributions, as prescribed by the model. 
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Q 2 


X 


y 


W 


e f f 


R 


F 2 


Ostat 


Osyst 


Stot 


(GeV 2 ) 






(GeV) 






(%) 


(%) 


(%) 


0.35 


0.0000061 


0.640 


240. 


0.384 


0.10 


0.397 


5.5 


14.5 


15.5 


0.50 


0.0000086 


0.640 


240. 


0.473 


0.13 


0.494 


3.6 


11.1 


11.7 


0.65 


0.0000112 


0.640 


240. 


0.539 


0.16 


0.568 


3.6 


10.4 


11.0 


0.65 


0.0000164 


0.440 


199. 


0.536 


0.15 


0.547 


2.9 


8.8 


9.2 


0.85 


0.0000138 


0.682 


248. 


0.664 


0.19 


0.713 


5.0 


13.0 


13.9 


0.85 


0.000020 


0.470 


206. 


0.628 


0.19 


0.646 


2.6 


6.6 


7.1 


0.85 


0.000032 


0.294 


163. 


0.615 


0.18 


0.621 


2.4 


7.2 


7.6 


0.85 


0.000050 


0.188 


130. 


0.576 


0.18 


0.578 


2.8 


8.1 


8.6 


0.85 


0.000080 


0.118 


103. 


0.533 


0.17 


0.534 


3.1 


15.4 


15.7 


1.20 


0.000020 


0.664 


245. 


0.793 


0.23 


0.857 


3.5 


10.5 


11.0 


1.20 


0.000032 


0.415 


194. 


0.741 


0.22 


0.759 


2.6 


7.1 


7.6 


1.20 


0.000050 


0.266 


155. 


0.709 


0.22 


0.715 


2.2 


6.9 


7.2 


1.20 


0.000080 


0.166 


122. 


0.626 


0.21 


0.627 


2.1 


8.9 


9.2 


1.20 


0.000130 


0.102 


96. 


0.569 


0.21 


0.570 


2.4 


5.1 


5.6 


1.20 


0.000200 


0.066 


77. 


0.525 


0.21 


0.525 


2.4 


5.2 


5.8 


1.20 


0.000320 


0.042 


61. 


0.531 


0.21 


0.531 


2.4 


7.9 


8.3 


1.50 


0.000032 


0.519 


217. 


0.817 


0.25 


0.855 


3.2 


8.8 


9.4 


1.50 


0.000050 


0.332 


173. 


0.771 


0.25 


0.783 


2.7 


5.9 


6.4 


1.50 


0.000080 


0.208 


137. 


0.687 


0.24 


0.690 


2.4 


6.9 


7.3 


1.50 


0.000130 


0.128 


107. 


0.668 


0.23 


0.669 


2.5 


6.6 


7.1 


1.50 


0.000200 


0.083 


87. 


0.645 


0.23 


0.645 


2.5 


7.4 


7.8 


1.50 


0.000320 


0.052 


68. 


0.613 


0.23 


0.613 


2.5 


6.9 


7.4 


1.50 


0.000500 


0.033 


55. 


0.577 


0.23 


0.577 


2.5 


7.9 


8.3 


2.00 


0.000032 


0.692 


250. 


0.945 


0.29 


1.048 


5.0 


12.9 


13.9 


2.00 


0.000050 


0.443 


200. 


0.908 


0.28 


0.939 


3.2 


6.0 


6.8 


2.00 


0.000080 


0.277 


158. 


0.743 


0.27 


0.751 


2.9 


5.8 


6.5 


2.00 


0.000130 


0.170 


124. 


0.727 


0.26 


0.729 


2.7 


7.6 


8.0 


2.00 


0.000200 


0.111 


100. 


0.716 


0.26 


0.717 


2.8 


5.8 


6.4 


2.00 


0.000320 


0.069 


79. 


0.727 


0.25 


0.727 


2.8 


5.9 


6.6 


2.00 


0.000500 


0.044 


63. 


0.639 


0.25 


0.639 


2.9 


8.3 


8.8 


2.50 


0.000050 


0.554 


224. 


0.973 


0.30 


1.034 


4.3 


11.2 


12.0 


2.50 


0.000080 


0.346 


177. 


0.932 


0.29 


0.950 


3.3 


5.4 


6.3 


2.50 


0.000130 


0.213 


139. 


0.917 


0.28 


0.922 


2.9 


7.6 


8.1 


2.50 


0.000200 


0.138 


112. 


0.840 


0.27 


0.842 


2.8 


6.8 


7.4 


2.50 


0.000320 


0.086 


88. 


0.702 


0.27 


0.703 


3.1 


5.3 


6.1 


2.50 


0.000500 


0.055 


71. 


0.649 


0.26 


0.649 


3.2 


6.8 


7.6 


2.50 


0.000800 


0.035 


56. 


0.590 


0.26 


0.590 


3.3 


8.2 


8.9 


3.50 


0.000080 


0.484 


209. 


1.045 


0.32 


1.094 


4.1 


8.6 


9.5 


3.50 


0.000130 


0.298 


164. 


1.003 


0.31 


1.018 


3.5 


5.0 


6.1 


3.50 


0.000200 


0.194 


132. 


0.869 


0.30 


0.873 


3.3 


7.2 


7.9 


3.50 


0.000320 


0.121 


105. 


0.898 


0.29 


0.899 


3.2 


8.5 


9.1 


3.50 


0.000500 


0.077 


84. 


0.863 


0.28 


0.864 


3.4 


6.5 


7.3 


3.50 


0.000800 


0.048 


66. 


0.686 


0.28 


0.686 


3.6 


7.4 


8.2 


3.50 


0.001300 


0.030 


52. 


0.663 


0.27 


0.663 


3.6 


7.5 


8.3 



Table 1: Proton structure function F2 (x, Q 2 ) and effective virtual photon-proton cross section 
a^lp(W, Q 2 ), scaled by the kinematic factor k = Q 2 / (47r 2 a), with statistical (S s t a t ), systematic 
(Ssyst ) an d total (5tot ) fractional errors. The normalization uncertainty, not included in the 
systematic error, is 3%. The values of R used to calculate F2 are also given. 
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The following contributions to the systematic errors from the event selection were in- 
cluded: trigger and timing veto 0.5%; BDC efficiency 2%; vertex finding efficiency 
2%. For the region y < 0.05 the systematic error on the vertex finding efficiency was 
increased to 5%. 

For the radiative corrections an error of 2% is taken everywhere, except for the highest y 
point of each Q 2 bin and for all points with Q 2 < 0.65 GeV 2 , where it is increased to 5% 
for the electron method. This error is due to uncertainties in the hadronic corrections, 
in the cross section extrapolation towards Q 2 = 0, in the higher order corrections and 
the absence of soft photon exponentiation in the HERACLES Monte Carlo. These 
effects were studied using the program HECTOR [41]. 

The uncertainty due to photoproduction background was assumed to be 30% of the 
correction applied, i.e. smaller than 6%. This affects only the highest y bins at low Q 2 . 

An additional error of 3% was assigned to the measurements using the £ method to allow 
for the uncertainties in the hadronic final state simulation of the Monte Carlo programs. 
This error was determined by comparing the results of different reconstruction methods 
for the hadronic final state, using the combined information of calorimeter cells and 
tracks, or by using calorimeter cells only. 

The effect of reduced efficiency for detecting diffractive events, such as the exclusive 
channel ep — ► epp°, has been estimated using the Monte Carlo programs for diffractive 
processes discussed in Section 4. Cross section corrections of up to 6% are applied 
for the points at the highest y values, and half of the correction was added to the 
systematic error. The effect is largest at the highest y values where the decay products 
of the meson often escape detection in the CJC and FT, and hence no event vertex is 
found. 

The overall normalization uncertainty is 3% due to the uncertainty in the luminosity 
determination. 



These systematic uncertainties affect differently the Fi measurements made with different 
methods. In Fig. || the comparison of the measurements made with the E and with the 
£ method is shown. The statistical and systematic errors are added in quadrature. The 
agreement between the two data sets is very good. 

For the final result the values obtained with the E method are takenf], supplemented with 
X points at low y (0.03 < y < 0.12) where no points from the E method are available. The 
result is shown in Fig. and tabulated in Table |l|. Both F2 and cr^lp are given. The latter has 
been multiplied by Q 2 / (47t 2 a) to demonstrate directly the effect of R. The table also contains 
the statistical, systematic and total errors and the value of R used for the F2 calculation. 
A table delineating the many different correlated and uncorrelated error contributions is 
available on request to the HI collaboration. The measurements have a typical systematic 
error of 5-10%. Compared with the previous HI analysis [Q] the F2 measurement has been 
extended to lower Q 2 (from 1.5 GeV 2 to 0.35 GeV 2 ), and to lower x (from 3 • 10~ 5 to 6 • 10~ 6 ). 

5 This procedure is used except for the highest x point at Q 2 = 1.5 GeV 2 . For this point the total error 
calculated with the £ method is almost a factor of two better than that of the E method, and therefore the 
£ method is used. 
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7 Discussion of the Results 



In Fig. |^ the F2 data are compared with previous HI measurements P], with the fixed target 
measurements of E665 [42| and NMC [39|, and with the predictions of models for F2 at low 
x. In the region of overlap the results are in good agreement with our previous measurements 
and the total error has been reduced by a factor of 2 to 3. The data also show a smooth 
continuation from the fixed target measurements towards the low-x region at HERA. The 
rise of F2 with decreasing x is still clearly prominent for values of Q 2 > 2 GeV 2 but becomes 
less steep for smaller Q 2 values. 

In Fig. ||, values of A are shown from fits of the form F2 oc x~ x (~ W 2X ) at fixed Q 2 to the 
present HI F2 data. For each Q 2 bin with Q 2 > 0.85 GeV 2 and x < 0.1 the exponent A was 
determined taking into account the point to point systematic error correlations. Note that 
this is merely a convenient parameterization to show the change of the slope with various Q 2 
in a quantitative way. Due to the HERA kinematics the x region in which the data are fitted 
is different for the Q 2 values used. The result is given in Table Fig. || also includes the 
HI measurements of A reported in [||. The newly measured values are consistent with the 
previous low-Q 2 measurements, and are in the range ~ 0.1 — 0.2. 



Q 2 (GeV 2 ) 


A 


5^stat 


fi^syst 


0.85 


0.146 


0.033 


0.101 


1.2 


0.192 


0.018 


0.038 


1.5 


0.128 


0.019 


0.036 


2.0 


0.133 


0.022 


0.038 


2.5 


0.216 


0.020 


0.043 


3.5 


0.172 


0.021 


0.031 



Table 2: The values of the exponent X as a function of Q 2 , together with the statistical 
(SXstat) an d systematic (5X sys t) contribution to the error. 



Fig. ^ shows function of the 7*p center of mass energy W for fixed values of Q 2 . 

The HI data are compared with low energy measurements and with photoproduction data. 
The low energy measurements are at slightly different Q 2 values and were propagated to the 
values indicated on the figure using the ALLM parameterization p3j (see below). 



In Fig. 10 crtt p is shown as a function of Q 2 for W values above 60 GeV. The data 
are transformed to the W values given using the ALLM parameterization. The new results 
presented in this analysis help to close the gap between measurements of real and virtual 
photoproduction and provide further information on the transition between real and virtual 
photon interactions. 

Several phenomenological parameterizations based on models have been proposed to de- 
scribe the region at low Q 2 and the transition from photoproduction to DIS, often using 
ingredients both from Regge theory at low Q 2 , and from QCD when Q 2 is of the order of 
1 GeV 2 or larger (for a recent review, see e.g. ||). A smooth transition from the photopro- 
duction to the DIS regime is generally assumed, but there is still uncertainty as to how this 
transition takes place and as to the underlying dynamics. 
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Parameterizations motivated by Regge theory relate the structure function to Reggeon 
and pomeron exchange phenomena which successfully describe the slow rise of the total 
cross section with the center of mass system energy in hadron-hadron and real photon- 
proton interactions. The model of Donnachie and Landshoff || (DOLA) assumes a pomeron 
intercept of 1 + A ~ 1.08 for the energy dependence of the cross section, independent of 
the virtuality of the photon. Fig. [j] shows that the model gives, as expected from the A 
measurements as a function of Q 2 , a prediction of F2 which is too low for the region of 
Q 2 > 0.85 GeV 2 , but which approaches the data for the lowest Q 2 values. The DOLA 
prediction can also be interpreted as the contribution of soft pomeron exchange to the cross 
section at non-zero Q 2 . 



In the approach of Capella et al. |LCj (CKMT) it is assumed that screening corrections 
and multi-pomeron exchange contribute less with increasing Q 2 , leading to a Q 2 dependent 
power A. Hence, the structure function F% predicted by CKMT rises faster with decreasing 
x for increasing Q 2 than the DOLA calculations. Furthermore, the model assumes that 
this prescription accounts for non-perturbative contributions to F2 at a Qq scale of about 2 
GeV 2 , from where, at higher Q 2 , perturbative QCD (pQCD) evolution equations are applied 
to predict the Q 2 dependence of F%. As expected, Fig. [?] shows that the prediction^ is 
systematically above that from DOLA, but it is still below the data for all but the lowest 
values of Q 2 . 



The parameterization of Abramowicz et al. p3| (ALLM) assumes that the total j*p cross 
section consists of two contributions which distinguish Reggeon and pomeron exchange and 
the behaviour of the power A is assumed to vary with Q 2 in a logarithmic way, emulating 
pQCD evolution in the high Q 2 region. In Fig. ^ it is shown that the prediction is below the 
data for small Q 2 values, but agrees with the data for Q 2 > 2 GeV 2 . 

Contrary to these Regge inspired models, the model of Gliick, Reya and Vogt |30| (GRV) 
is defined completely within the parton picture. It is assumed that all parton distributions 
at very low Qq = 0.34 GeV 2 have a valence-like shape, i.e. vanish for x — > 0. Assuming 
further that F2 evolves towards large Q 2 values via parton radiation which is given by the 
leading twist QCD evolution equations, GRV predict that the structure function F2 should 
rise with decreasing x, even for low values of Q 2 ~ 1 GeV 2 . Fig. ^ shows that the GRV 
distributions describe the data well for Q 2 > 1 GeV 2 , but are systematically lower than the 
data for Q 2 < 1 GeV 2 . At low Q 2 values the QCD evolution is over a range of Q 2 which is 
too small, and the parton distributions become dominated by the valence behaviour at the 
starting scale Qq. 



The model of Badelek and Kwiecinski [38] (BK) combines the concepts of Vector Meson 



Dominance (VMD) with dynamical parton models such as that of GRV. It has per force 
a smooth transition from pQCD to the real photon limit, which coincides with the region 
measured here. The BK model predicts the F2 values very well as shown in Fig. |?]. In Fig. || 
BK is compared with the total j*p cross section measurement and also shows good agreement. 

A different approach to the low Q 2 behaviour of F2 in the transition region between 
photoproduction and DIS has been presented in |ll[ by Adel, Barreiro and Yndurain (ABY). 
It assumes that perturbative QCD evolution is applicable to the lowest values of Q 2 . In 
addition to the flat or so called "soft" behaviour of F2 with decreasing x for Q 2 ~ 1 GeV 2 , 



6 The CKMT curves shown in Fig. (?] were calculated using the parameters given in |Io) without QCD 
evolution. Therefore the prediction is not shown for the highest Q 2 value. 
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which is manifest in for example the GRV model, a "harder" contribution is introduced to 



prevent F2 decreasing with decreasing x for Q 2 values below 1 GeV 2 . In [11] the prescription 
for this contribution is singular ~ x~ Xs , with X s = 0.48 independent of Q 2 . Furthermore the 
strong coupling constant is assumed to become independent of Q 2 for values below roughly 
1 GeV 2 , that is a s "saturates". The result of a fit of the ABY approach to data, reported 
11 1, is compared with the new F2 results presented here in Fig. [?]. There is good agreement 



m 



at low x, but also some possible disagreement in the higher (x, Q 2 ) region shown in the figure. 
Fig. |?| also shows that in this prescription, at low Q 2 , the rise of F2 with decreasing x occurs 
at lower x compared to the other approaches, namely for x < 1CP 4 . 

In Fig. pi], BK and ALLM are compared with the measurements of cr^l^. To ensure 
consistency, for the BK prediction the value of R has been taken from [37]. The ALLM 
prediction as calculated here is for R = and therefore it is not shown for the largest W. 
For the other values of W the effect of R is a few percent. The measurements suggest that 
in the current ALLM parameterization the transition towards photoproduction behaviour 
occurs at too large a Q 2 , in contrast to the BK approach. However BK predict a photopro- 



duction cross section which is larger than the measurements, as shown in the figure. Fig. 1C 
also shows the HI fit based on BK, as discussed in Section 6. The photoproduction cross 
section measurements at HERA were used in the fit, which leads to a stronger turn-over 
to the photoproduction regime compared to the original BK model. This fit gives a good 
phenomenological description of all data shown in the figure. 

In summary, it turns out that the region 0.1 < Q 2 < 1 GeV 2 spans the kinematic range 
in which Regge or VMD inspired models describe the data at low Q 2 , and models based 
on pQCD account well for the higher Q 2 domain. Although in general several of these 
models show the correct qualitative behaviour of the data, none of them gives at present a 
completely satisfactory description of the data. Future studies will have to show if this can be 
remedied by adjusting the parameterizations whilst simultaneously preserving their internal 
consistency. The data in the low (x,Q 2 ) region will help to discriminate between different 
theoretical approaches to low-Q 2 dynamics. 

Finally a parameterization is given of HI data in the region 0.35 < Q 2 < 5000 GeV 2 
and y > 0.01. Starting from the double asymptotic expression for F2 as developed in ||44fl , a 
QCD inspired parameterization with only two parameters was reported in f45| . Here, F2 is 
parameterized as 



NfX -wvHm (8) 



with T = ln[ln((Q 2 + Qg)/A 2 // )/ln(Qg/A 2 // )] and N f = 571/^77^/324, where n f is the 
number of quark flavours, taken to be equal to four. The parameter A e yy resembles the QCD 



mass scale. The function T has been slightly changed with respect to the original proposal |4J 
so as to achieve a smooth behaviour over the whole Q 2 range. The parameter 7 is defined as 
^12/(11 - 2n f /3). Fitting this expression to the HI data gives Ql = 373 ± 24 (stat.) ± 42 
(syst.) MeV 2 and A e// = 247 ± 11 (stat.) ±19 (syst.) MeV, with a x 2 /ndf = 156/231 (full 
errors). Athough eqn. |8| does not satisfy the constraint, imposed by current conservation, 
that F2 should vanish at Q 2 = 0, it constitutes a compact Q CD-inspired parameterization of 
all HI measurements of F2. 
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8 Summary 



A measurement has been presented of the proton structure function F2(x,Q 2 ) and virtual 
photon-proton cross section a^lp(x,y,Q 2 ) in deep inelastic electron-proton scattering at low 
Q 2 with data taken in the 1995 HERA running period. These are the first measurements 
made with the upgraded backward calorimeter and drift chamber of the HI detector. 

The measurements presented are obtained using two different methods to reconstruct the 
inclusive scattering kinematics, allowing both a powerful internal cross check of the data 
and the measurement in a large kinematic region. The data cover the region of Q 2 between 
0.35 GeV 2 and 3.5 GeV 2 and with Bjorken-x values down to 6 • 1CP 6 . The measurements 
show a smooth transition from the fixed target high-x data to the HERA low-x data. 

The distinct rise of the structure function with decreasing x in the low-x region, which 
is very prominent for Q 2 > 2 GeV 2 , diminishes at lower Q 2 values. When taken together 
with the data from fixed target experiments, the rise observed for the smallest Q 2 values 
approaches that expected in Regge and VMD interpretations. 

The data have been compared with different models which aim to describe the whole Q 2 
region. Several of these models predict the correct qualitative behaviour observed in the data 
but presently do not agree with the data throughout the full kinematic range. The data 
access the transition region from DIS to photoproduction and provide powerful constraints 
on the development of further low-Q 2 phenomenology. 
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Figure 1: Experimental (points) and Monte Carlo (solid lines) distributions of a) the energy 
in the electron tagger, and b ) the energy of the (fake ) electron candidate in the SPACAL, for 
tagged photoproduction events. 




Figure 2: Experimental (points) and Monte Carlo (solid histograms) distributions of a) 
the energy of the scattered electron and b ) the polar angle of the scattered electron for DIS 
events. The Monte Carlo curves are the sum of DIS and photoproduction events and the 
photoproduction events alone. 
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Figure 3: Experimental (points) and Monte Carlo (solid histograms) distributions ofY>(Ei — 
P z ,i) measured in the calorimeter for a) all DIS event candidates and b) DIS event candidates 
with y e > 0.55. The Monte Carlo curves are the sum of DIS and photoproduction events and 
the photoproduction events alone. 




Figure 4: Experimental (points) and Monte Carlo (solid histograms) distributions for y% > 
0.05 of the ratios of a) the y values measured with the S and E method yz/Ue, and b) the 
transverse momentum of the hadronic system and the electron Pth/Pt,e- The Monte Carlo 
curves are the sum of DIS and photoproduction events and the photoproduction events alone. 
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Figure 5: Comparison of the proton structure function F%{x, Q 2 ) as a function of x at various 
values of Q 2 (in GeV 2 ) measured with the E method (full points) and with the £ method 
(open points). The errors represent the statistical and systematic errors added in quadrature. 
A global normalization uncertainty of 3% is not included. 
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Figure 6: Measurement of the proton structure function F2(x,Q 2 ) as a function of x at 
various values of Q 2 (in GeV 2 ). The inner error bars are the statistical errors, the outer 
error bars represent the statistical and systematic errors added in quadrature. A global 
normalization uncertainty of 3% is not included. 
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Figure 7: Measurement of the proton structure function F2(x, Q 2 ) in the low Q 2 region by HI 
(full points), together with previously published results from HI (open circles), E665 (open 
triangles), NMC (open squares). The Q 2 values are given in GeV 2 . Various predictions for F2 
are compared with the data: the model of Donnachie and Landshoff (dashed line), the model 
of Capella et al. (dotted line/small ), the model of Badelek and Kwiecinski (dashed-dotted 
line), the model of Gliick, Reya and Vogt (full line) and the model of Adel et al. (dotted 
line/large). Global normalization uncertainties are not included. 
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Figure 8: Variation of the exponent A from fits of the form ~ x~ x at fixed Q 2 values and 
x < 0.1. Full symbols are the data from this analysis; open symbols are the data from ||/. 
The inner errors are statistical, and the full errors represent the statistical and systematic 
errors added in quadrature. 
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Figure 9: Measurement of the total virtual photon-proton cross section (J^*p 3>S <% function 
ofW at various values of Q 2 (in GeV 2 ). The cross sections are multiplied with the factors 
indicated in the figure (numbers in brackets). The errors represent the statistical and sys- 
tematic errors added in quadrature. Full symbols are data from this analysis; open symbols 
are previously published data from HI (circles), E665 (triangles) and NMC (squares). Global 
normalization uncertainties are not included. The curves represent the DOLA (dashed line, 
only shown for Q 2 = 0), ALLM (dotted line) and BK (dashed-dotted) parameterizations. 
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Figure 10: Measurement of the virtual photon-proton cross section Cy*p clS Q, function of Q 2 at 
various values ofW (in GeV). The cross sections for consecutive W values are multiplied with 
the factors indicated in the figure (numbers in brackets). The errors represent the statistical 
and systematic errors added in quadrature. Full symbols are data from this analysis; open 
symbols are HI data from The photoproduction points (cross: W = 210 GeV, diamond: 

W = 170 GeV) are from Global normalization uncertainties are not included in the 

errors shown. The curves represent the ALLM (dotted line) and BK (dashed-dotted line) 
parameterizations and the HI fit based on BK (dashed line). 
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